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Abstract 

This article deals with higher order Caputo fractional variational 
problems with the presence of delay in the state variables and their in- 
teger higher order derivatives. 



1 Introduction 

Lagrangian theories involving higher-order derivatives appear of great interest 
as an imbedding for field theories with fields of subcanonical dimension, e.g. 
the Heisenberg nonlinear spinor theory, for which local interactions are less 
singular than in the canonical case. 

The calculus of variation has a long history of communications with other 
fields of mathematics such as geometry and differential equations, and with 
physics. Recently the calculus of variations has found applications in economics 
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and some branches of engineering such as electrical engineering. Optimal con- 
trol which is a rapidly expanded field can be regarded as a part of the calculus 
of variations. 

Recently, the fractional calculus which is as old as the classical calculus has 
become a candidate in solving problems of complex systems which appear in 
various branches of science PII^ISllill ^ lTlfn i fn]. 

Several authors found interesting results when they used the fractional cal- 
culus in control theory (531 [51] . 

Experimentally, the use of delay together with the fractional calculus may 
give better results. 

Optimal control problems with time delay in calculus of variations were dis- 
cussed in |25) . Variational optimal control problems within fractional deriva- 
tives were considered in [2 7) . Fractional variational problems in the presence 
of delay were studied in [SI [TU| . 

The aim of this paper is to deal with optimal control fractional variational 
problems in the presence of delay in the state variable and its higher order 
derivatives 

The structure of the paper is as follows: 

In section 2 the necessary definitions in the fractional calculus used in 
this manuscript are reviewed. In section 3 the unconstrained fractional Euler- 
Lagrange equations with delay are discussed . The fractional control problem 
is presented in Section 5. 

2 Basic Definitions 

We present in this section some basic definitions related to fractional deriva- 
tives. 

The Left Riemann-Liouville Fractional Integral and The Right Riemann-Liouville 
Fractional Integral are defined respectively by 



where a > 0, n — 1 < a < n. Here and in the following T{a) represents the 
Gamma function. 

The Left Riemann-Liouville Fractional Derivative is defined by 




(1) 





The Right Riemann-Liouville Fractional Derivative is defined by 



1 / rl \ ^ 

WW^r>^(-s) / (-r-'-'/M-i., (4) 
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The fractional derivative of a constant takes the form 



(6) 



and the fractional derivative of a power of t has the following form 

T(p - a + 1) 

for /3 > -l,a > 0. 

The Caputo's fractional derivatives are defined as follows: 

The Left Caputo Fractional Derivative 
and 

The Right Caputo Fractional Derivative 

where a represents the order of the derivative such that n — 1 < a < n. By 
definition the Caputo fractional derivative of a constant is zero. 

The Riemann-Liouville fractional derivatives and Caputo fractional deriva- 
tives are connected with each other by the following relations: 

fc=o a + ij 

^D^it) = Dtf{t) Xi ViT-^^^^ (10) 
t^r(fc-a + l) 

In [T], a formula for the fractional integration by parts on the whole interval 
[a, h] was given by the following lemma 

Lemma 2.1. Let a > 0, p,q > 1, and ^ + ^^<l + a(p^l and q ^ 1 in 
the case when + ^ — 1 + a ) 

(a) If (f G Lp{a, h) and tp £ Lq{a, b), then 

^{t){al^m)dt = (\{t){I^^){t)dt (11) 
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(b) Ifg e I^{Lp) and f G aI"{L,), then 

\{t){aD^f){t)dt^ f f{t){D^g){t)dt, (12) 



where aI"{Lp) := {f ■ f = al"g,g e Lp{a,b)} and 
IbiLp) :={f:f = I^g,geLp{a,b)}. 

In [S] and [TD] , other formulas for the fractional integration by parts on the 
subintervals [a, r] and [r, b] were given by the following lemmas 

Lemma 2.2. Let a > 0, p,q > 1, r ^ {a, b) and j^ + ^<l + a(py^l and 



P 9 

(a)I{ (fi G Lp{a, b) and tp G Lq(a, b), then 



q ^ 1 m the case when j^ + ^~l + a) 



J a 

and thus if g G Iri^p) f ^ aI^{Lq), then 

g{t){ aD^m)dt = f f{t){D'^g){t)dt (14) 

J a 

(h)If (p G Lp{a^ b) and ip G Lq{a, b), then 

J r J r 

ilj{t){( ip{s){s-t)°'~^ds)dt (15) 



r(a) J a 

and hence if g G I^{Lp) and f G aI°'{Lq), then 



g{t){ aD^f){t)dt = / f{t){Dtg){t)dt 

{aD^f){t){l {Dtg){s){s~tr-^ds)dt. (16) 
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That is 



r(«) Ja 

g{t){aD^f){t)dt= I f{t){D^g){t)dt 
1 



r(a) Ja 



b 



jm-i iD^g)is)is-tr-^ds)dt. (17) 
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Lemma 2.3. Let a. > 0, p,q > 1, r G (a, b) and ^ + ^ < 1 + a ( p ^ 1 and 

q ^ \ in the case when - -\--~l-\-a). 

(a) I{ (fi e Lp{a, b) and ip S Lq{a, b), then 

(\{t)m){t)dt= fmiri^^mdt (is) 

J r J r 

and thus if g S rI"{Lp) and f £ I^{Lq), then 

f g{t){D'Sf){t)dt = f f{t){rD'^g){t)dt (19) 

J r J r 

(b) If (f € Lp{a^ h) and ijj € Lq(a^ h), then 

f\{t){bi''m)dt = f m{i:vmdt 

J a J a 

^f{ajf ~ sr-'ds)dt (20) 

and hence if g G aI"{Lp) and f G I^{Lq), then 

f g{t){DU){t)dt = f f{t){aD"g){t)dt 

- [{DUmi j\aD-g){s){t - sr-'ds)dt. (21) 

f g{t){Dtm)dt = f mi^D^gmdt 

~f{ajf fit)rD''ij\aD^9){s){t - sr-Hs)dt. (22) 

3 The Unconstrained Caputo Fractional Vari- 
ation with delay 

Before we consider the fractional control problem, let us consider the following 
fractional variational problem with delay 
Minimize 

J{y)= f L[t, ^D"^y{t), ^D^^y{t), '^D^-y{t), 

J a 

^£>f^2/(t), ^£>f=2/(t),..., ''D^,-y{t),y{t),y'{t),...,y^''Ht), 



That is 
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y(i-r),y'(i-T),...,2/«(i-r)]dt, 

such that 



(23) 



fc - 1 < amax < k, amax = maxjai, /3j}i<i<„.i<i<m, 

y(')(&) = Q, / = 0,1,2, ...,fc- 1, y{t)^(b{t) te [a-T,a], 

T>0,T<b-a, Q!,,/3j e R V z = 1,2, Vj = 1, 2, m, (24) 

where c; are constant, (f>{t) is a smooth function and L is a function with con- 
tinuous first and second partial derivatives with respect to all of its arguments. 

If the above variational problem (j23p has a minimum at yait) and ri(t) S R 
is an admissible function such that rjit) = in the interval [a — r, a] then the 
function 

at) = J(2/o + ev) (25) 
where e £ R admits a minimum at e = 0. Hence 

«6 n m 

/ 5,+iL(t) a^i?«'77(t) + ^ a„+,+iL(i)^i?f^?7(t)+ 
i=i j=i 

k k 
J2 dm+n+p+2m7J^P\t) + ^ d,n+n+k+3mTJ^P\t - r)]dt = 0. (26) 

where diL is the partial derivative of L with respect to its argument. 
On using the connection formulas © and (flUl) . reads 

" 2=1 j=l 

k k 

J2 dm+n+p+2L{t)ii^P\t) + ^ a™+„+fc+3£(t)r7'''ni - r)]dt = 0. (27) 

p=0 p=0 

Now if one splits the integral, makes the change of variables for i — r and 
uses the fact that 77 = in [a — r, a] , ([27]) becomes 

/ a.+iL(0 aZ?"'r;(0 + ^ dn+,+iL{t)Dl' ^{t) + 

i=i j=i 

fc fc 
^ 5„+„+p+2£(t)r/(P) (t) + dm+n+k+-3L{t + ^)^'''^ 

p=0 p=0 
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/I) n 1 1 L 

--^ 1=1 j=i 

k 

J2 dm+n+p+2mr,^P^ {t)]dt = (28) 

By using the integration by parts formulas in the mentioned above Lemma 
12. 1[ Lemma 12.21 and Lemma 12.31 and the usual integration by parts formula , 
one obtains the following 



i=i 

dp _ ^„dp 



p=0 p=0 
n 1 nb 

YT-^^b-ri / {D'^,^d,+,L){s){s - tr-^ds)Ut)dt+ 

/I) n m 
[Y,D^^{d,+,L){t)+J2 b-rD^^{dn+,+iL){t)- 
— ^ ^=l j=l 



ds)]-q{t)dt+ 



dP 



^^"^^'d^ id-^+n+p+2L){tMt)dt+ 



p=0 
k p—1 



M 

EE(-1)V id^^n^p+2L){tW-^-\t)t-^ + 
p=l g=0 

E E(-i) V (^"+"+fc+3^)(^ + ^)^^"'"' 

EE(-l) V id^+n+p+2m)^P~''~\t)\l_^ = (29) 

In equation (|29|) if one chooses 77 such that 77(a) = and 77 = on [6 — t, 6], 
one gets 

n m 
E^"-.(9.+li)W +E ai?''^(a„+, + lL)(t) + 

i=l J=l 
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(5™+„+p+2i)(i)+E(-l)'rf^ {dm+n+k+r.+,L){t + r)- 

p=0 p=0 



n 1 j-b 



= 0. (30) 



In equation (|29p if one chooses 77 such that 77*^'^ (6) = and = on [a, 6 — t], 
one gets 

n m 

Y,D^^{d,+ ,L){t)+Y. b-rD^^{dn+j + lL){t)- 
i=l j=l 
rb-7 



1 nO — T 



E^"^^'d^ (a™+„+p+2i)W = 0. (31) 

p=0 



Now since both integrals in (|29l) are now zero, one gets 

EE(-i)'d^ (a™+„+p+2L)(i)r/^-'^-^(i)ir^+ 

p = l g=0 

k p-1 

p=l q=0 

EE(-l)''d^^ (a™+„+p+2L)(t)r?^-^-^(t)lt. =0. (32) 

p=l q=0 

Thus one can state the foUowing theorem 

Theorem 3.1. Let J(y) be a functional of the form i(23)) defined on a set 
of continuous functions y{t) which have continuous Caputo fractional left order 
derivatives of orders Oii and right derivative of order /3j in [a, b] and satisfy 
the conditions in {21 Let L : [a - t, b] x ^'"+"+2^+2 ^ ^ have continuous 
first and second partial derivatives with respect to all of its arguments. Then 
the necessary conditions that J{y) possesses a minimum at y{x) are the Euler- 
Lagrange equations 



3=1 ^b 



8 



" 1 riT 



for a < t < b — T, 



^ dP dL 



1 pb — T or 

5^ -— - / ( ,D^^{-—^){s){t ~ = (34) 

,^1 ^ IPjJ Ja ( 



'd^Dl^yit) 

for b — T < t < b, and the transversality conditions 

k p—1 

dti 



p=l q=a 



p=l q=0 ^ ^ ' 

p— 1 q—O ^ ' 

for any admissible function rj satisfying ri{t) = t E [a — t, a], 
rj^^^b) =0, / = 0, l,2,...,fc- 1. 

Theorem 13.11 can be generalized as follows 

Theorem 3.2. Consider the functional of the form 

J(yi,y2,...,yd) = / L[t, ^Z?"^yiW, ^i?"^yi(i), 

J a 

'^D'^^y.it), ^D"^y2{t),..., ^i^""y2 W, 
^D'^^y.it), ^D'^^y^it),..., ^D'^^y^it), 
^D^yi(t), ^D^,^y,{t),..., ''D^.-y.it), 

yi ) , ) , . . . , 2/f ^ (i ) , y2 ) , 2/2 ) , • ■ • , ) , • ■ • , 
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yd{t),y'd{t), ■■■,y''d\t),yi{t - T),y[{t - T),...,y[''\t - r), 

y2{t - r), y'.,(t - r), yf (t^r), y^t - r), y'^it - r), - T)]di, (36) 

jdeGned on sets of continuous functions yi{x), i = 1,2,..., d that have left 
Caputo fractional derivatives of order ai £ M., i ~ 1, 2, n and right Caputo 
fractional derivatives of order Pj G M, j = l,2,...m in the interval [a,b] and 
satisfy the conditions 

yf\h) = Cii, Z = 0, 1, k, yi(t) = (j)i{t) i = 1,2..., d, t e[a-T,a] 

a < b, T > 0, T <b - a , (37) 

where k-1 < Umax < k, amax = max{ai, /3j}i<i<,i,i<j<m, cu 's are constant 
and (pi 's are smooth functions and 

L : [a — T, 6] X Rrf(™+"+2fc+2) R jg a function with continuous first and second 
partial derivatives with respect to all of its arguments. For yi{x), i — 1,2, ...,d, 
satisfying if3^) to be a minimum of i|36)) . it is necessary that 



p=o 9yr (i) p=o dyr'(t-T) 

)(s)(s-i)"--ids) = (38) 



for a < t < b — T, z — 1,2, d 

for b — T < t < b, z — 1,2, d and the transversality conditions 

k p-1 
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p=i9=o '^^ dy'/'{t) 

for any admissible vector function rj = (771, 772, ?7d) satisfying 
rjit) = (0,0,...,0) t G [a-r,a], - (0,0, ...,0), d = 0, 1, 2, fc - 1. 

4 The Fractional Optimal Control Problem 

Find the optimal control variable u{t) which minimizes the performance index 
Jiy,u)= / F[t,u{t), ^D'^^yit), ^ D''^y{t), ^D''-y{t), 



^Z?f^y(t), ^i^f^y(t),..., ^DP-y{t),y{t),y'{t),...,y(''Ht), 

y{t~T),y'{t-T),...J^\t-T)]dt, (41) 
subject to the constraint 

G[i,«(<), ^i?"Hj(i), ^Z?"=2/(t),..., ^i5""2/(t), 
^Df^y(t), ^i?f^2;(<),..., ''D^,-y{t),y{t),y'{t),...J^\t), 

y{t~r),y'{t-T),...J^\t-T)]=Q (42) 

such that 

y(')(5) = Q, 0, l,2,...,fc- 1, j/(i) = 0(i) te[a-T,a], a<b,T>0 

T <b-a, a, e R, i = 1,2, ...,n, /Sj G M, j = 1,2, ...,to, (43) 
where c; are constant and F and G are functions [a — r, fe] x ]U"+"i+2fe+3 _^ jg^ 
with continuous first and second partial derivatives with respect to all of their 
arguments arguments. 

To find the optimal control, one defines a modified performance index as 

J{y,u) = I F + X{t)G dt, (44) 



where A is a Lagrange multiplier or an adjoint variable. Using the conditions 
([59)1 and (PD|) in Theorem 13. 21 the following necessary equations for opti- 
mal control are found: Euler-Lagrange equations 

i=l ^ \ / j^--^ 
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h h 



W + = 0, 



' ' ' du{t) 
for a < t < b — T, 

^ CI 77' ^ fi/^ 

^ (P , dF , A. 



du{ty ' ' 'au(t) 

for 6 — T < t < 6, and the transversality conditions 

^ jq pi TP 

p=l q=0 ^ ^ ^ 
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p=l g=0 ^ ^ ^ 

p=l 9=0 " ^ ' 

where 77 is any admissible function satisfying r]{t) = t & [a — T,a], 
T]^^\b) =0, ; = 0, 1,2, ...,A;-1. 

5 Conclusion 

In this manuscript we have developed a fractional control problem in the pres- 
ence of both left and right Caputo fractional derivatives of any order and delay 
in the state variables and their derivatives. The results were applied in or- 
der to find the necessary conditions for the optimal control problems. When 
ai 1, = l,...,n and I3j — )• 1, = l,...,m the classical problem is 
recovered. 
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